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Abstract 

Polar wander may occur on Triton and Pluto because of volatile migration. Triton, 
with its low obliquity, can theoretically sublimate volatiles (mostly nitrogen) at the rate of 
~10 14 kg yr" 1 from the equatorial regions and deposit them at the poles. Assuming Triton to 
be rigid on the sublimation timescale, after ~10 5 years the polar caps would become large 
enough to cancel the rotational flattening, with a total mass equivalent to a global layer 
-120-250 m in depth. At this point the pole wanders about the tidal bulge axis, which is the 
line joining Triton and Neptune. Rotation about the bulge axis might be expected to disturb 
the leading side/trailing side cratering statistics. Because no such disturbance is observed, it 
may be that Triton’s mantle viscosity is too high but its surface volatile inventory is too low 
to permit wander. On the other hand, its mantle viscosity might be low, so that any 
uncompensated cap load might be expected to wander toward the tidal bulge axis. In this 
case, the axis of wander passes through the equator from the leading side to the trailing side; 
rotation about this wander axis would not disturb the cratering statistics. Low-viscosity 
polar wander may explain the bright southern hemisphere: this is the pole which is 
wandering toward the equator. In any case the “permanent” polar caps may be 
geologically very young. Polar wander may possibly take place on Pluto, due to its 
obliquity oscillations and perihelion-pole geometry. However, Pluto is probably not 
experiencing any wander at present. The Sun has been shining strongly on the poles over 
the last half of the obliquity cycle, so that volatiles should migrate to the equator, stabilizing 
the planet against wander. Spacecraft missions to Triton and Pluto which measure the 
dynamical flattening could give information about the accumulation of volatiles at the poles. 
Such information is best obtained by measuring gravity and topography from orbiters, as 
was done for Mars with the highly successful Mars Global Surveyor. 
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1. INTRODUCTION 

Polar wander may occur on Triton and Pluto due to the migration of volatiles. 
Because of Triton’s low obliquity, the Sun shines strongest on the areas near the equator. 
Insolation falling on Triton’s equatorial regions is theoretically strong enough to sublimate 
the equivalent of a global layer ~ 120-250 m thick on a timescale of ~10 5 years and deposit 
the material at the poles (see Figure 1). This is enough to cancel the flattening due to 
rotation, so that if the satellite is rigid, Triton overbalances and the pole wanders. 

The pole can also wander due to the formation of polar caps if Triton’s effective 
viscosity is low. Consider the Earth. The Earth’s pole wanders due to density anomalies 
caused by tectonics. These anomalies are not nearly large enough to overcome the Earth’s 
equatorial bulge, as they would need to be in the rigid body case. Instead the bulge relaxes 
fast enough so that it migrates with respect to a coordinate system fixed in the surface; the 
excess matter associated with the density anomalies eventually ends up at the equator, which 
is the position of maximum dynamical stability. The mantle relaxation time is about 10 3 
years, which is short compared to the 10 8 year tectonic timescale; hence the Earth rapidly 
adjusts to the tectonic situation, and the pole wanders. It could be the same for Triton and 
for Pluto, with the wander being caused by the pile-up of volatiles, rather than tectonics 

The situation for Pluto is more complicated than for Triton because Pluto’s 
obliquity is such that the long-term insolation at the equator is nearly the same as at the 
poles. If polar wander does occur, it may be associated with the obliquity oscillations and 
the perihelion-precession geometry. The corresponding figures for Pluto seem likely to be 
a global layer ~90-180 m thick sublimated on the obliquity timescale of ~10 6 years. 
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2. TRITON 

Triton orbits Neptune. Triton’s orbital plane is tilted with respect to Neptune’s 
equator by inclination I = 156.8°, so that the motion is retrograde. Triton is in synchronous 
rotation about Neptune, with the satellite’s equatorial plane being nearly coincident with its 
orbital plane. Moreover, Neptune’s equator is inclined its own orbital plane about the Sun 
by e N = 29.6°, where £ N is Neptune’s obliquity. As a result of this geometry, the subsolar 
point on Triton is always confined to a band between latitudes ±53° (Harris, 1984; Yelle et 
al., 1995), so that the equatorial regions receive more insolation than the poles. Thus, in 
addition to seasonal frost deposition (Hansen and Paige, 1992; Spencer and Moore, 1992), 
there will also be a secular trend to sublimate volatiles (principally nitrogen) in the equatorial 
regions and condense them at the poles (Brown and Kirk, 1994). This secular trend could 
have consequences for polar wander, the gravitational field, resurfacing rates, and cratering 
statistics. 

2.1 Insolation and sublimation on Triton 

Assuming Neptune’s orbit about the Sun is circular, the long-term latitudinal second 
degree (/= 2) A F component of the insolation incident on Triton can be written as 


fa Vf 


AF= F t 


\ U N J 


16 


|(^sin 2 £ r -£)p 2 (sin0) 


( 1 ) 


(e.g., Rubincam, 1994), where F E = 1371 W m' 2 is the solar constant at the Earth distance of 
a E = 1 AU from the Sun, a N = 30.1 AU is Neptune’s distance from the Sun (AU = 
Astronomical Unit), § is latitude on Triton, £j. is Triton’s obliquity relative to the Sun, and 


4 



4 


Rubincam 6/18/02 


5 


P 2 (sin ())) = (3 sin 2 <|> - l)/2 

is the second degree Legendre polynomial. The obliquity &j. constantly changes because 
Triton’s orbit precesses on Neptune’s equatorial plane with a period of 600 to 700 years 
(e.g., Harris, 1984; Jacobson, 1990). The relationship is 

cos e T = cos I cose N + sin / sin e N cos Q\t -t 0 ) (2) 

where Q' is the precession angular speed, t is time, and t 0 is some reference time. Plugging 
this equation into (1) and averaging over time yields 

<AF> = F e («£/a N ) 2 [(5/16) P 2 (sin <j>) (2 - 6 cos 2 / cos 2 £p - 3 sin 2 / sin 2 e,. )/8] (3) 

which, for the values of / and ^ given above, yields 

<AF> = - 0. 1 25 P 2 (sin <(>) W m' 2 (4) 

so that there is excess insolation in the equatorial regions, as stated above. For a typical area 
of 2k Rj (where R T = 1352 km is the radius of Triton), an insolation value of -0.125 W m' 
2 \ an albedo A = 0.8, and a latent heat of sublimation of 2.5 x 10 5 J kg' 1 (Spencer et al., 

1997, p. 454), the rate of sublimation S is approximately 

(1 - 0.8) (2) (3.14) (1.352 x 10 6 m) 2 (0.125 W m' 2 ) (3.16 x 10 7 s yr' 1 )/(2.5 x 10 5 J kg' 1 ). 


or an amount 
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5 = 3.6 x 10 13 kg yr (5) 

can be sublimated from the equatorial regions for redeposition at the poles. This is much 
greater than the ~10 9 -10 10 kg yr' 1 produced by the plumes observed so far on Triton (Kirk 
et al., 1995). Also, (5) works out to be approximately 36 km 3 yr' 1 , which is much greater 
than the resurfacing rate conservatively estimated at 0.01 km yr' 1 by Stem and McKinnon 
(2000). Hence the rate of volatile movement can theoretically be quite high. 

2.2 Mass distribution and gravitational Held 

The next step is to find how this sublimation might affect the mass distribution on 
Triton. The mass distribution is most conveniently expressed in terms of Triton’s 
gravitational field, which can be written in the usual spherical harmonic form 

V(r,Q,X) = S P ^(sin^)(QmCOsmA + 5 to sinmA) (6) 

r £=oV T J m=0 

(e.g., Kaula, 1966; Caputo, 1967), where V(r,<|>,A,) is the gravitational potential, r is the radial 
distance from the center of mass, <j> is latitude, X is east longitude, and G = 6.673 x 10' 11 m 3 
kg' 1 s' 2 is the universal constant of gravitation, while A/ T = 2. 14 x 10 22 kg is the mass of 
Triton. The P fm (sin <| >) are the associated Legendre polynomials, while the C tm and S (m are 

unnormalized spherical harmonic coefficients. The convention among geodesists is to 
speak of the zonal C m coefficients in terms of J e , where simply J ( = -C e0 . 

Triton feels centrifugal forces from rotation, with the rotational potential being 
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f (o\ r 2 \ 

3 J 


[P 2 (sin <p) + 2] 


(7) 


where P 2 (sin (|)) = P 20 (sin <|>). This potential tends to flatten the satellite, making Triton 
somewhat pumpkin-shaped. Triton’s responds to this potential with J 2 tol , the second degree 
spherical harmonic coefficient due to the rotational hydrostatic flattening (e.g., Kaula, 1968, 
pp. 68-73): 


_ C0 2 R}. 

= 3 GM t 




-1 . 


( 8 ) 


In the above equation co r = 2 tt/( 5.877 days) is the angular speed of Triton’s rotation about 
its axis, while 


At = CjIMjRj 2 (9) 

is the central condensation factor, with C T being its polar moment of inertia. Here At 
measures how much differentiation has taken place; it should not be confused with the A 
which is used for longtiude in (6). At present C T is unknown, but substituting At = 0.4, the 
value for a homogeneous sphere, yields by (8) a maximum value for J 2 rot of 

J 2 rot * 1.325 x 10 4 , A,. = 0.4 . (10) 
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Rotation is not the only contributor to the gravitational field; the bulge raised by 
Neptune’s tidal attraction will also give rise to C tm and S lm coefficients. Neptune’s tidal 

potential is 


V M , = GM ? - P,(c osT), (11) 

a 

where T* is the angle measured from the line joining Triton and Neptune, M H = 1.024 x 10 26 
kg is Neptune’s mass, and a = 3.548 x 10 8 m is Triton’s distance from Neptune. The tidal 
potential tends to stretch Triton into a shape something like an American football, with the 
long end pointing to Neptune. It is interesting to note that the ratio of the coefficients of the 
Legendre polynomial in (7) and (1 1) is (tides/rotation) = 3GM n ldyfci'~ 3. The tidal 
potential is thus 3 times stronger than the rotational potential, so that most of Triton’s 
flattening is due to tides and not rotation. In other words, Triton has a distinctly triaxial 
shape, being more of a football than a pumpkin, with the axis of the football in the equatorial 
plane. 

The C fm and S lm tidal coefficients can be found with the help of the addition theorem 
(e.g., Kaula, 1968, p. 67): 

P,(cos¥) = £ <2 ~ y i : — P,.(sinfe)P,.(sm^)cos[m(A - A„)] . (12) 

~ 0 (£ + m)\ 

Here S y is the Kronecker delta, defined such that 8 y = 1 for i = j and is 0 otherwise. In the 
above equation <J) N and are the respective latitude and east longitude of the sub-Neptune 
point. Apart from libration, Neptune will be on Triton’s equator, and by convention 
longitude is measured from the sub-Neptune point. Libration will be ignored, so that <j) N = 0 
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and A^ = 0 for the purposes considered here. Thus substituting (12) into (11) yields by 
analogy with (8) 


^ 2, tide 


Icv’Mj 


5 



15 

4 




(13) 


for the equilibrium tides, which gives 

J 2Me = 1.986 x 10' 4 , ^ = 0.4 (14) 

for the J 2 component of the tidal bulge raised on a homogeneous Triton by Neptune. The 
flattening from the tidal potential is thus larger than the rotational flattening. The total 
hydrostatic coefficient will be 

“^2, hydro = ^2,rot ^2, tide ~ 3.311 X 10 , Aj = 0.4 (15) 

It may be that J 2 rot is smaller than given by (15) because Triton is differentiated, as 
would be expected if it was captured into orbit by slamming into previously existing moons 
and was greatly heated by that and the tidal circularization of its orbit (Goldreich et al., 
1989). If Triton is highly differentiated, as seems likely from the collision scenario, then if 
the rocky core has a constant density p r and the icy mantle a constant density p i; the 
condensation factor (9) can be shown to be (Rubincam, 2000) 

fa — - X = — (Pr ~ £j ) Pj_ 

T M T Rl 5 [Pt (p r - Pi) 213 p T _ 
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where p T = 2064 kg m' 3 is Triton’s mean density (e.g., McKinnon et al., 1995). Water ice 
with density p, = 1000 kg m' 3 and anhydrous rock with density p r = 3500 kg m' 3 give ^ = 
0.3 1 1. On the other hand, for hydrated rock with density p r = 2800 kg m' 3 (McKinnon et 
al., 1997, p. 302), A r could only get as low as 0.339. Assuming a worst case value of ^ = 
0.300 gives 


and 


J 2mt - 0.645 x lO 4 , 

Ap«0.3 

(16) 

■A, tide ~ 0.966x 10\ 

Aj- ~ 0.3 



.hydro ^2. rot + y 2 , tide =i.6iix 10 


,-4 


Ay « 0.3. 


(17) 


For principal axis rotation C 21 = S 2l = 0. Also, S 22 = 0 by (12) because 7^ = 0. 
Thus the only other second degree coefficient will be C 22 , given by 


= RjM n 

' 22Mde ~ 4 a 3 M T 


i |5 15, 

1 + A, 

2 4 7 


-1 


(18) 


which gives 


C 22 , tide = 0.993 x IQ 4 , Aj = 0.4 (19) 


for a homogeneous Triton, and 
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C 22tide = 0.483 x 10- 4 , \ = 0.3 (20) 

for a highly differentiated Triton. The numbers computed here are summarized in Table 1. 

2.3 Rigid body (high interior viscosity) case 

Uncompensated polar volatile deposits will act to decrease J 2 because they add mass 
to the poles, making the satellite dynamically rounder. This section calculates the time 
necessary for a rigid Triton to overcome the rotational flattening. A more rigorous but less 
physically obvious derivation is relegated to the Appendix. 

Assuming that the ice deposit has the shape of a spherical cap centered on the pole 
with a mass A/ cap , angular extent \j/, and has a uniform thickness, the contribution to J 2 will 
be 


^2, cap = “^cap ( C0S V + ^S 2 VV ( 2M t) ( 2 1 ) 

(Rubincam, 1992, p. 2630). Assuming a somewhat arbitrary angular size of \j/ = 45° makes 
the above expression become 


h cap * “^cap/(2M T ) . (22) 

with smaller cap sizes giving bigger contributions to J 2 . 

If the effective viscosity of Triton’s mantle is high, so that Triton can be considered 
rigid on the sublimation timescale, then the polar deposits must overcome the rotational 
flattening for polar wander to occur. This happens when 
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4™ +4, „ = 0 . (23) 

The polar caps cancel the rotational flattening and dynamically only the tidal bulge is left; 
Triton is then free to rotate about the Neptune-Triton line, i.e., the pole is free to wander 
about the axis of the tidal bulge, as shown in Figure 2. Any uncompensated mass centered 
on the pole will be in unstable equilibrium, so that any perturbation will cause Triton to 
rotate to try and place the excess mass on the equator. Thus the time t s required for 
sublimation to overcome the rotational flattening will, by (5), (10), (22), and (23), be given 
by 


*s = S = 2MjJ 2 ro , / S ~ 1.6 x 10 5 yr , ^ = 0.4 (24) 

for a homogeneous Triton, where the amount of volatiles involved is 

M cap « 6xl0 ,8 kg, ^ = 0.4 (25) 

which is equivalent to a global layer 250 m deep, assuming a density of 1000 kg m' 3 . For a 
highly differentiated Triton the numbers are 


t s = MJ S = 2MjJ 2 rot / S * 0.8 x 10 s yr , X, = 0.3 (26) 


and 


M cap «3xl0' 8 kg, ^ = 0.3. (27) 
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The above figures are the equivalent of a global layer only 120 m deep and a timescale of 
just ~8 x 10 4 yr. On the other hand, there may be some nonhydrostatic (e.g., tectonic or 
impact) component of unknown size to overcome, so that t s and M cap might be larger than 
given above. These results are summarized in Table 1. 

2.4 Low interior viscosity case 

The above values for t s and M cap assumed Triton to be a rigid body. In reality, the 
surface will relax under the load and J 2 will decay according to a characteristic time. For 
a homogeneous viscous planet the / = 2 Darwin relaxation time t D toward isostatic 

compensation is (Darwin, 1879): 


t D = 19iy(2p T £A) ( 28 ) 

where r) T is Triton’s effective viscosity, p T = 2064 kg m' 3 is Triton’s mean density, and g T = 
0.78 m s' 2 is the surface gravitational acceleration (e.g., McKinnon et al., 1995). The 
viscosity r) T is not known. For comparison, an Earth-like mantle viscosity of r| T = 10 21 Pa s 
yields 


t D = 1 .38 x 10 5 yr, r| T = 10 21 Pa s (29) 

and the relaxation time is on the order of the sublimation time t s . 

Isostatically compensated mass anomalies “float” on the substrate like icebergs and 
have only a small effect on the gravity field and moments of inertia. Thus the sinking of the 
ice load toward isostacy inhibits polar wander. On the other hand, a low interior viscosity 
means that the rotational and tidal bulges can relax quickly, so that the pole will wander 
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southward in order to place whatever uncompensated ice load there is on the equator. This 
situation is more complicated than the rigid body case and will not be treated in detail here, 
except to make two remarks. 

The first remark is that in the low viscosity case the ice load need not overcome the 
rotational flattening for the pole to wander, as it did in the rigid body case, so that the 
uncompensated ice cap need not be as thick. This is analogous to the situation on the Earth: 
plate tectonics causes the pole to wander, even though the mass associated with tectonics is 
nowhere near large enough to overcome the Earth’s huge rotational flattening. On the Earth 
values of C 2m and S 2m having typical magnitudes of 10' 6 , while J 2 = 10' 3 . But the Earth’s 
rotational bulge relaxes fast enough to allow the pole to shift, thanks to a mantle viscosity 
which gives a t D = ~10 3 yrs, which is short relative to the 10 8 yr tectonic timescale (e.g., 
Stacey, 1992, pp. 394-400; Lambeck, 1980, pp. 342-346). Likewise, ice ages on the Earth 
create ice caps which cause a certain amount of polar wander, even though these caps are 
nowhere near the size of the equatorial bulge (e.g., Bills and James, 1997). 

The second remark is that in the low viscosity case, Triton’s pole might be expected 
to wander in the direction shown in Figure 3, so that Triton attempts to place the extra mass 
not just on the equator, but on the axis of the tidal bulge, assuming that there are no other 
nonhydrostatic loads. The reason for putting it there is that this is the position of maximum 
dynamical stability, just as the extra mass on the Moon’s equator is at the sub-Earth sub- 
anti-Earth axis. The sense of rotation is 90° away from the rigid body case. If there are 
such nonhydrostatic loads, as might be created by convection in Triton’s interior (analogous 
to what happens on Earth), then it is not clear what the polar wander path is. It might be a 
random walk, as it is on Earth; or if the ice load overwhelms the other sources, it may still be 
predominately as shown in Figure 3. Low-viscosity polar wander may be the reason for the 
large observed south polar cap, as discussed later (section 4). 
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3. PLUTO 

Pluto can be treated in a way analogous to Triton. As with Triton, the migration of 
volatiles to Pluto’s poles can have consequences for resurfacing rates, polar wander, the 
gravitational field, and cratering statistics. 

3.1 Pluto insolation 

For Pluto the expression analogous to (1) is 

<AF>= ^fe) (M! sin!£f 4h (sinw (30) 

where a p = 40 AU is Pluto’s distance from the Sun and e P is Pluto’s obliquity. Pluto’s 
obliquity oscillates with time according to the equation 

e P = 115.5°+ 11.8° sin x (31) 

where % = 19.5° + 130.2° t, with time t being measured in millions of years from the 
present time (Dobrovolskis and Harris, 1983; Dobrovolskis et al., 1997). Thus Pluto’s 
obliquity oscillates between -103.7° and -127.3° with a 2.8 x 10 6 period. For obliquities 
less than -125° the poles receive more insolation than the equator; see Ward (1974) and 
(30) above. Therefore Pluto spends nearly all of its time with the poles receiving more 
sunlight than the equator (e.g., Dobrovolskis et al., 1997; Spencer et al., 1997). 

A naive interpretation of this situation would be that over the long term, volatiles flee 
to the equator, making the planet if anything flatter instead of rounder, so that there is no 
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polar wander on Pluto. However, it may be that Pluto’s climate dynamics do not decree that 
volatile migration goes to zero at the obliquity where the polar insolation exactly equals the 
equatorial insolation (Rubincam, 2000). Perhaps instead of occuring at e p = 125°; volatile 
migration reverses course nearer its present average obliquity of 1 15.5°. It might even be 
that Pluto has altered its average obliquity through the phenomenon of climate friction so 
that it is now in a state of no net transport of volatiles over geologic time; however, it is by 
no means clear at present whether climate friction can have either the magnitude or the 
direction required to perform this feat (Rubincam, 2000). 

If Pluto’s average obliquity is near a state of zero net transport, then volatiles may 
periodically shuttle into and out of the polar regions as the obliquity oscillates. Volatile 
movement may also depend on the perihelion-spin pole geometry. Pluto’s orbit about the 
Sun is highly eccentric, with an orbital eccentricity of 0.25. Pluto’s spin pole precesses on 
its orbital plane with a period of 2.8 x 10 6 years (Dobrovolskis et al., 1997). As Pluto 
precesses, it exposes first its equator and then the poles to the most sunlight at perihelion. 
The question becomes whether during times of polar deposition volatiles accumulate to the 
point where the pole wanders. 

In the absence of a detailed climate model, a typical insolation deficit will be 
assumed to be a factor of 10 lower than for Triton: 

AF = - 0.0 1 P 2 (sin <|>) W nf 2 (32) 

giving, in a calculation analogous to (4), 


S = 2.1 x 10 12 kg yr 1 , 


(33) 


using R p = 1150 km as the radius of Pluto. 
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3.2 Rigid body (high viscosity case) 

The rotational and tidal flattenings can be computed as was done for Triton. Just as 
Neptune raises a tidal bulge on Triton, so Charon raises a tidal bulge on Pluto. Using the 
values M p = 1.288 x 10 22 kg, M c = 1.891 x 10 21 kg, © p = 271/(6.387 days), and a = 19,662 
km in (8) and (13) yields 


7 2rot = 1.146xl0-\ V = 0-4 (34) 

and 

J 2Me = 0.220 x 10' 4 , ^ = 0.4 (35) 


leading to 


A**. = A„ + Aw. = 1 - 3 « 6 * 10' 4 , K = 0.4 (36) 

It will be seen that the hydrostatic tidal flattening relative to rotational flattening is much 
smaller for Pluto than for Triton, but is still significant. These numbers lead to 

t s = 2MpJ Xlol / S ~ 1.6 x 10 6 yr , A* = 0.4 (37) 

by analogy to (16), and giving 

M cap = 3xlO ,8 kg, Ap = 0.4 (38) 
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for a rigid Pluto. This is equivalent to a global layer 180 m deep. Since t s is about half the 
obliquity oscillation period and the precession period, polar wander seems barely possible if 
Pluto is rigid, assuming the supply of volatiles is great enough for the planet to overbalance. 
For a highly differentiated Pluto, the corresponding figures are 

J 2rot = 0.558 x 10' 4 , Xp = 0.3 (39) 

and 

J ’ 2 ude = o. 107 x 10- 4 , \,= 0-3 (40) 

leading to 

^2, hydro = 0.665X 10" 4 , ^ = 0.3 (41) 

t s = 2MpJ 2 wt / S = 8 x 10 5 yr , /ip = 0 - 3 (42) 

and giving 

M cap = 2xl0 18 kg, \ = 0.3. (43) 


This last value is equivalent to global layer whose thickness is 85 m. 

The C 22 tide term can be calculated for Pluto by using the analog of (18). The result 
is 
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C 22itide = 0.110 Xp = 0.4 (44) 

C 22jlide = 0.054 = 0.3 . (45) 

These results are summarized in Table 1. 

3.3 Low viscosity case 

For Pluto p p = 2022 kg m' 3 , g p = 0.65 m s' 2 , so that the Plutonian equivalent of (28) 
is 

h> = 19rip/(2pp gpR p ) (46) 

This gives 

t D « 2 x 10 5 yr, ri P = 10 21 Pa s (47) 

where 10 21 Pa s is an Earth-like mantle viscosity. Unlike the case for Triton, t D is shorter 
than the assumed value for t s if Pluto has an Earth-like viscosity. Marcialis (1985) has 
argued that Pluto in fact has a very low interior viscosity, so that the planet relaxes very 
quickly, yielding a topographically bland planet. 

The direction of Pluto’s polar wander will be analogous to that of Triton. In the 
rigid body case it will be as shown in Figure 2; while for a low viscosity it will presumably 
be as shown in Figure 3, assuming no other nonhydrostatic loads. 
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4. DISCUSSION 

For an Earth-like mantle viscosity of 10 21 Pa s, Triton’s relaxation time t D is about 
the same as the sublimation time t s ; while for Pluto t D is shorter than t s . The actual mantle 
viscosities of both Triton and Pluto are unknown. McKinnon et al. (1997, p. 316) use a 
viscosity estimate of t) = 139 exp(A/7) Pa s for Pluto’s water ice mantle, where A = 6830 K 
and T is measured in kelvins. This equation gives rj > 10 21 Pa s for the outer -150 km of 
the planet for a particular thermal model, but obviously the viscosity is very sensitive to the 
chosen value of A. 

Triton’s total planetary inventory of nitrogen has been estimated at a global surface 
layer equivalent of -1 km (Brown and Kirk, 1994; Yelle et al., 1995, p. 1068). This is close 
to the -120-250 m necessary for rigid body wander to occur, so that a substantial fraction 
of all its nitrogen would have to be on the surface. Also, Triton’s cratering statistics argue 
against rigid body wander having occurred over the age of the current surface. Triton’s 
leading side (apex of motion) has more craters than its trailing side (Croft et al., 1995). If 
Triton’s interior viscosity is high compared to the sublimation timescale, then the pole 
would be expected to wander in the sense shown in Figure 2, which would tend to smear out 
the leading side/trailing side asymmetry and produce a symmetrical ring whose axis is the 
Triton-Neptune line, assuming the pole wanders continuously in one direction. Hence the 
current data seem consistent with Triton having a high mantle viscosity but with a low 
surface volatile supply, the supply being too small to overcome the flattening. However, it 
may be that Triton rocks back and forth about the bulge axis instead of rolling 
continuously; in such a case the smearing would be less. 

On the other hand, the cratering statistics may be compatible with polar wander 
taking place on a low viscosity interior. In this situation, the pole would be expected to 
wander in the sense shown in Figure 3 if the caps are the only uncompensated loads, so that 
the axis of wander pierces the leading side and trailing side on the equator. Rotation about 
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this axis will not disturb the leading side/trailing side asymmetry. Another argument for a 
low viscosity is that Triton could be very active geologically (Stem and McKinnon, 2000), 
which may mean that the interior is hot and its viscosity relatively low. Certainly if Triton 
has a subsurface ocean the effective viscosity will be extremely low. 

It must be stated that only 40% of Triton’s surface has been seen clearly enough to 
yield useful cratering statistics (Stem and McKinnon, 2000). The other half of the satellite 
may yield something interesting. 

If polar wander does occur, then ice in the caps is brought closer to the equator, 
where sunlight can once again sublime the material for transport back to the poles. On 
Triton deposition is expected to be more or less continuous, so that the pole will always be 
wandering. This will have implications for the thickness of the ice caps and for resurfacing 
the satellite. 

If Triton can be considered rigid on the sublimation timescale, then for two ice caps 
each with angular size of 45°, the ice will cover the topography to an average depth no 
greater than -850 m. The corresponding number for Pluto is -615 m. If the mantle 
viscosity of either body is low, then presumably the depths are even smaller. In any case 
there will be a limit to the depths of the ice caps on the order of hundreds of meters. 

Brown and Kirk (1994) propose the viscous spreading of the caps to explain why 
Triton has no thick ice caps concentrated at the poles, but polar wander seems a more likely 
explanation: the deposits move south and are sublimated before they can form thick sheets. 
The speed at which the pole wanders is not computed here, but the timescale would be 
expected to be short, so that the “permanent” polar caps would be geologically very young. 

Polar wander may be the explanation for the anomalously bright southern 
hemisphere: perhaps it is this hemisphere which is wandering toward the equator, bringing 
the bright surface deposits with it, while the northern hemisphere tilts back out of view. 
Another possibility is that the bright feature is a seasonal polar cap (Yelle et al., 1995). 

(“North” and “south” here follow the convention of referring to these directions 
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on Neptune and not Triton. But under the normal assumption that the positive spin vector 
pierces the northern hemisphere, the bright hemisphere is actually Triton’s northern 
hemisphere. Recall Triton’s motion is retrograde.) 

A spacecraft mission to Triton might give information about the moon’s stiffness 
and the migration of volatiles. Suppose polar wander occurs continuously on Triton; then 
for a rigid interior, the moment of inertia about its rotation axis would always be near the 
limit of stability, and Triton’s gravitational field would be dramatically different from what 
is expected at hydrostatic equilibrium due to rotation and tides. For instance, if the 
spacecraft measured a J 2 much lower than the / 2 h . given by (17) for a highly 
differentiated Triton, then this would be dramatic evidence for polar wander. 

A spacecraft measurement of Pluto’s J 2 would be less clear-cut. Pluto’s obliquity 
has just passed through a period of “low” obliquity, “low” in this case meaning values 
closer to the minimum numerical value of 104° rather than the high of 127°, as can be seen 
from (31). Maximum insolation on the poles occurs at 90° obliquity and decreases at 
higher values. So the Sun has been shining as much as it ever does on Pluto’s poles for the 
last ~1.5 x 10 6 years or so. Hence volatiles might be expected to sublimate and the poles 
accumluate at the equator, stabilizing the planet against wander and leading to a greater than 
average value of J 2 at the present time. But excess mass from any other source, such as 
tectonic activity, would be expected to end up at the equator and increase J 2 , just as is the 
case on Earth. Hence there would be a problem in determining from spacecraft how much 
of the excess mass is due to volatile migration and how is much due to other sources. 

Whether polar wander occurs on Pluto due to volatile migration is problematic. For 
the simple climate model given by (30), volatile accumulation at the poles takes place only 
when the obliquity is near its extreme value of 127.3°. The value for S given by (33) is an 
extremely crude estimate, and assumes (30) is not quite the right model. A better climate 
model is needed for Pluto to see whether the wander is actually possible. 
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Should one or the other of the gravitational fields be unexpectedly quiet, with only 
hydrostatic components present indicative of a fluid body, then the moment of inertia could 
be found from inverting (18) as well as the sum of (8) and (13). This yields 
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using the ^22, tide coefficient for Triton, with corresponding equations for Pluto, giving 
estimates on the degree of differentiation inside Triton and/or Pluto. For instance, for the 
Earth the observed value for J 2 is 1.0826 x 10 3 (e.g., Lerch et al., 1979; Stacey, 1992, p. 
136), which by (45) leads to Ag = 0.332, a figure which is basically indistinguishable from 
that obtained from precession (e.g., Stacey, 1992, p. 33). The reason the numbers are so 
close is that the Earth’s rotational flattening dwarfs all other sources of the gravitational 
field. On the other hand, for Mars J 2 = 1.959 x 10' 3 (Smith et al., 1993), yielding A^ = 
0.375. The observed value from precession is A^, = 0.365 (Folkner et al., 1997); which is 
smaller than the inferred value. The difference is due to a large nonhydrostatic component 
from Tharsis and other sources. It may be that similar estimates can be made for Triton and 
Pluto without being too badly corrupted by other sources. 

A spacecraft orbiting Triton and Pluto would give the best determinations of J 2 , C 22 , 
and S 22 , as opposed to a flyby mission. A laser altimeter aboard an orbiter would also 
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provide topographic and gravitational field data, just as in the case of Mars Global Surveyor 
(Zuber et al., 1998; Smith et al., 1999). 

Hammel et al. (2001) propose a Neptune mission which might make repeated 
flybys of Triton. While not as useful as a Triton orbiter, such a mission would still provide 
valuable information about the satellite’s gravity and topography. 

A Pluto orbiter would be more involved. Unlike the proposed Neptune mission, 
which uses Neptune’s atmosphere to slow down, aerobraking cannot be used at Pluto 
because its atmospheric envelope is too tenuous. The spacecraft’s speed would have to be 
killed with rocket thrust, making the spacecraft heavier and more expensive. 

5. CONCLUSIONS 

Polar wander may occur on Triton due to the insolation-driven migration of 
volatiles. Triton’s apparent large south polar cap might be caused by its wander toward the 
sub-Neptune point. This style of polar wander would imply that Triton’s mantle has a 
relatively low viscosity. A low viscosity would likely agree with the leading side-trailing 
side cratering statistics and probably agree with the observed internal activity. Polar wander 
on Pluto due to volatile migration is more problematic, but still possible. Spacecraft 
measurements of Triton’s gravitational field could give evidence for or against volatile- 
induced polar wander on that body. Spacecraft could possibly measure the degree of 
differentiation inside Triton and Pluto. 

Appendix 

This appendix gives a slightly more rigorous derivation of (24) via the power 
balance equation. The insolation and deposition of nitrogen at the surface will be assumed 
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to be governed by the simple equation (e.g., Hansen and Paige, 1992; Spencer and Moore, 
1992; Yelle et al., 1995): 


(1 - A)F = e ss c SB r - L-j- + k-^- (Al) 

at dz 

where T is Triton’s surface temperature, A is its albedo, e ss is its emissivity, and k its thermal 
conductivity Also, F is the insolation, o is the mass per unit area of solid nitrogen on the 
ground (measured in kg m' 2 ), z is vertical direction (positive downwards), t is time, while L = 
2.53 x 10 5 J kg' 1 is the latent heat of sublimation for nitrogen and a SB = 5.67 x 10' 8 Wm' 2 K' 
4 is the Stefan-Boltzmann constant. 

It will be assumed here that Triton’s nitrogen atmosphere is always on the vapor 
pressure curve, and that solid nitrogen covers the satellite’s surface completely; in this case 
Tis a constant across the satellite (Yelle et al, 1995). It will also be assumed that A and k 
are constant across the surface as well. The insolation F and surface mass distribution a 
can be expressed in terms of spherical harmonics: 
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For the / = 0, m = 0 term, (Al) becomes 
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For all other terms, (Al) is 


(1 ~A)F? m = 



(A5) 


0-A)F f s m = 


dt 


(A6) 


Temperature drops out because T is taken to be constant across the surface; it is therefore 
found only in the / = 0, m = 0 term given by (A4). 

Only the /= 2, m = 0 term in (A5) affects J v In (A2) the insolation coefficient is 
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(Rubincam, 1994). Thus integrating (A5) gives 



5(1 -A)F e 
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where by (2) and (3) in the main text 


<3 sin 2 &,. - 2> = (1 - 3 cos 2 / cos 2 £ N - 3 sin 2 / sin 2 £ N /2) . 


The gravitational flattening due to the surface mass distribution can be shown to be (e.g., 
Kaula, 1968, p. 66): 
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A-- C M - (5)^' (A9) 

where p T = 2064.3 kg m' 3 is the average density of Triton. 

To overcome the rotational flattening 

^2, rot /z.vo! = 0 • (A 10) 

so that the time t s required for sublimation to overcome the rotational flattening will, by 
(A8)-(A10), be 
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This equation gives sublimation times almost identical to those given by (24) and (26). 
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Table 1 


X 

J 2, rot 

^ 2, tides 

*^2, hydro 

r 

'-'22, tide 

d 

h 

Triton 

(10- 4 ) 

(10- 4 ) 

(10- 4 ) 

(10- 4 ) 

(m) 

(10 5 yr) 

0.3 

0.645 

0.966 

1.611 

0.483 

120 

0.8 

0.4 

1.325 

1.986 

3.311 

0.993 

250 

1.6 

Pluto 

0.3 

0.558 

0.107 

0.665 

0.054 

85 

8? 

0.4 

1.146 

0.220 

1.366 

0.110 

180 

16? 


Parameters for Triton and Pluto. X = C/MR 2 , where C is the moment of inertia 
about the rotation axis, M is mass, and R is radius, d is the equivalent depth of the global 
layer of solid N 2 needed to make an ice cap which cancels the total hydrostatic flattening J 2 
hydro on a rigid object, where J 2 hydro = J 2 , m , + J 2 , ,i de • h is the time needed to sublimate the 
equivalent global layer and condense it at the poles in order to overcome J 2 hydro . 
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Figure 1 

The Sun shining on the equatorial regions of Triton could cause volatiles 
(principally nitrogen) to migrate to the poles. 

Figure 2 

Path of polar wander on a rigid body. The ice caps alias the rotational flattening, 
leaving only the tidal bulge. The axis of wander is then the axis of the tidal bulge. Rotation 
about this axis would be expected to smear out leading side-trailing side cratering statistics. 

Figure 3 

Expected path of polar wander on body whose mantle has a low effective viscosity. 
The pole will wander so as to attempt to place the excess mass of the caps on the axis of the 
tidal bulge. The polar wander path is 90° away from that shown in Figure 2, with the axis 
piercing the leading side and trailing side on the equator. Rotation about this axis does not 
disturb the leading side-trailing side cratering statistics. 
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Polar Wander on Triton and Pluto 
Due to Volatile Migration 

by 

David Parry Rubincam 


POPULAR SUMMARY 

Polar wander may occur on Triton and Pluto because of volatile migration. Triton, 
with its low obliquity (obliquity is the tilt of the axis), can theoretically sublimate volatiles 
(mostly nitrogen) at the rate of ~10' 4 kg yr' ! from the equatorial regions and deposit them at 
the poles. Assuming Triton to be rigid on the sublimation timescale, after ~10 5 years the 
polar caps would become large enough to cancel the rotational flattening, with a total mass 
equivalent to a global layer -120-250 m in depth. At this point the moon overbalances and 
the pole wanders about the tidal bulge axis, which is the line joining Triton and Neptune. 
Rotation about the bulge axis might be expected to disturb the leading side/trailing side 
cratering statistics; the leading side bashes into more meteorites and asteroids. Because no 
such disturbance is observed, it may be that Triton is too stiff but its surface volatile 
inventory too low to permit wander. On the other hand, its stiffness might be low, so that 
any size cap might be expected to wander toward the tidal bulge axis. In this case, the axis 
of wander passes through the equator from the leading side to the trailing side; rotation 
about this wander axis would not disturb the cratering statistics. Polar wander may explain 
the bright southern hemisphere: this is the pole which is wandering toward the equator, 
while the other hemisphere tilts out of view. In any case the “permanent” polar caps may 
be geologically very young. Polar wander may also take place on Pluto, due to its obliquity 
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oscillations and perihelion-pole geometry. However, Pluto is probably not experiencing any 
wander at present. The Sun has been shining strongly on the poles over the last million 
years, so that volatiles migrate to the equator, stabilizing the planet against wander. 
Spacecraft missions to Triton and Pluto which measure how flat the objects are could give 
information about the accumulation of volatiles at the poles. Such information is best 
obtained by measuring gravity and topography from orbiters, as was done for Mars with the 
highly successful Mars Global Surveyor. 
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